Morphological Symmetry Breaking During Epitaxial Growth 

at Grazing Incidence 
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It is shown that, in submonolayer growth at off-normal incidence, even much less than one percent 
of transfer from the condensation energy of the deposited atoms into adatom motion is sufficient 
to induce a net adatom current from the illuminated edge of a two-dimensional island to the other 
edges, thereby breaking the island symmetry. Such a symmetry breaking phenomenon is most 
pronounced for deposition at grazing incidence. Comparison between our theoretical predictions 
and existing experimental results confirms the general validity of the model. 



PACS numbers: 68.55.Jk, 68.35. Bs, 68.70+w 

Because of their application potentials in future elec- 
tronic devices, a great deal of effort is being devoted to 
developing novel methods for fabrication of organized 
low-dimensional structures, such as ordered arrays of 
quantum wires or quantum dots Q]. Epitaxial growth 
and its inverse process, atom removal by sputtering or 
etching, are two of the most promising approaches for 
mass production of controlled nanostructures on various 
substrates. Typically, ordered structures are obtained in 
heteroepitaxial systems, and the long-range elastic field 
associated with the lattice mismatch between the two sys- 
tems plays an essential role in leading to self-organized 
growth 0], However, in recent studies of Cu(100) ho- 
moepitaxy, an intriguing phenomenon about island for- 
mation and ordering has been discovered, i.e., the average 
island symmetry varies with the incident angle of depo- 
sition, leading to the formation of elongated islands, or 
ripples, of twofold symmetry |3| . In contrast, earlier stud- 
ies of Cu(100) homoepitaxial growth at normal incidence 
resulted in only square-shaped islands of fourfold sym- 
metry [0JjJ. Such a symmetry breaking phenomenon is 
already present even in the submonolayer growth regime, 
and is more dramatic at grazing incidence P|. For atom 
removal by sputtering under ion bombardment at off- 
normal incidence, formation of coherent ripples has also 
been observed on different substrates . 

In earlier attempts to understand these symmetry 
breaking phenomena, some qualitative suggestions have 
been proposed, all of which relying on atoms climbing 
down from steps as an essential atomic process [|]||. In 
this Letter, through a detailed study of a simple model, 
we offer an alternative and quantitative interpretation 
of the widely observed incidence geometry induced sym- 
metry breaking. Our model is based on an experimen- 
tally widely-invoked concept, namely, transient mobility 
of deposited atoms on various surfaces [|8| — |37T[] . Through 
studying shape evolution of monolayer-high islands on 



an fcc(lOO) surface, we demonstrate that even much less 
than one percent of transfer of the condensation energy 
from the deposited atoms into adatom motion is sufficient 
to induce a net adatom current transferring adatoms 
from the illuminated island edge to neighboring edges. 
As a consequence, the symmetry of the growing islands 
changes from initially square shape to rectangular shape 
elongated perpendicularly to the incident direction. Such 
a symmetry breaking phenomenon is most pronounced 
for deposition at grazing incidence. A comparison be- 
tween theory and experiment confirms the general valid- 
ity of the model. We also make several specific predic- 
tions that can be tested in future experiments. 

We use Fig. 1 to schematically show the growth pro- 
cess of a monolayer-high island on an fcc(lOO) surface 
during deposition at off-normal incidence. We denote 
the island width and length by W = m x a and L = m y a, 
respectively, where a is the surface lattice constant. We 
name the four island edges as x + , x~ , y + , and y~ edges. 
Atoms are deposited onto the surface at an angle a with 
respect to the surface normal, with a deposition flux F. 
We consider the three most important kinetic processes 
for adatom motion, i.e., the surface diffusion, the island 
edge diffusion, and the island corner crossing, with rates 
qi = Vi exp (— Vi/kT), where Vi and (i = s,e,c) are 
the corresponding barriers and attempt frequencies. In 
general, one has V c = V e + AV, with AT^ > because 
an adatom has to lower its coordination in crossing an 
island corner |r^ , [l3|| . 

In the case of deposition at off-normal incidence, the ef- 
fective flux for depositing atoms on the surface is F cos a. 
The adatoms diffuse on the surface and attach to an ex- 
isting island with a flux f s — F ™ s - , where N is the 
island density. Following classical nucleation theory [O] , 
we have 
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where 6 is the coverage. Because of the isotropic nature 
of adatom diffusion on an fcc(fOO) surface, the fluxes for 
adatom attachment at the x (x + or x~) and y (y + or 
y~) island edges are given by f x = 2 { m 7+m v ) fs and h = 
tjt — — r/s, respectively. A remarkable feature for the 
off-normal deposition is that it directly deposits atoms on 
the island edge facing the incident beam (the illuminated 
edge), with a flux fa = Fm y a 2 sin a, valid in the low 
coverage limit. 




FIG. 1. Schematic illustration of island growth on an 
fcc(lOO) surface at off-normal incidence (see text for nota- 
tions) . 

It is easy to conjecture on the island shape evolution 
if without the contribution of the condensation energy. 
Because of the existence of the direct landing flux fa, 
the illuminated edge (the y + edge in Fig. 1) has more 
arriving atoms than the other three edges. In the low 
temperature regime where island-corner crossing is not 
frequent enough to establish equilibrium adatom distri- 
bution along the four edges, the islands would elongate 
parallel to the incident direction. In contrast, at higher 
temperatures, frequent edge diffusion and corner cross- 
ing can take place, leading to adatom transfer from one 
edge to the other edges. This yields a uniform adatom 
distribution on the four island edges. Therefore, without 
the contribution of condensation energy, one would only 
find compact islands of either square shape or rectangular 
shape elongated parallel to the incident direction. 

Now, let us consider the effect of the condensation en- 
ergy transfer to adatom diffusion. As stated above, the 
illuminated island edge has more landing atoms than the 
other three edges. The existence of the transient mobility 
|p|-|rT| implies that those adatoms directly deposited at 
the illuminated edge will have higher mobility to cross the 
two island corners bounding this edge. Such anisotropic 
island corner crossing induces a net current transferring 
adatoms from the illuminated edge to its neighboring 
edges, and may therefore change the symmetry of the 
initially square-shaped islands. As a result, the islands 
can elongate along the direction perpendicular to the in- 
cident direction, which is just the opposite to the expec- 



tation given above without consideration of condensation 
energy transfer. 

To quantify the above picture, next we derive a set of 
equations describing the island evolution in the presence 
of condensation energy transfer. In a unit time, the total 
condensation energy given up by the atoms directly de- 
posited at the illuminated edge amounts to faUo, where 
Uq is the condensation energy for each atom. During the 
same time, the number of atoms landing on the illumi- 
nated edge is fa + f y . Because only a portion of the con- 
densation energy is transferred into diffusional motion, 
we write the average energy gain of an adatom on the 
illuminated edge as AE = (3faUo/(fa + fy)- Here (3 is a 
parameter describing the transfer efficiency from the con- 
densation energy of the incident atom to its diffusional 
motion, with (3=1 corresponding to complete transfer 
and (3 = to zero transfer. It should be noted that, 
besides its being finite based on the experimental evi- 
dences for transient mobility [p|-pT|, very little is known 
about this parameter (3. Nevertheless, it is known to pos- 
sess the following qualitative features, based on simple 
physical considerations. First, (3 is weakly temperature- 
dependent for a given system, because the condensation 
energy is much larger than the thermal energy at typical 
growth temperatures. Second, (3 is also weakly depen- 
dent on the incident angle, because the primary angular 
dependence of the problem has already been incorporated 
into the flux expression fa. Thirdly, (3 is expected to 
be strongly system-dependent, larger for a larger mass 
mismatch between the incident atom and the substrate 
atom. This last point is transparent within the simple 
picture of elastic collision between the incident atom and 
a surface atom, and is still qualitatively correct even if 
the frictional forces due to vibrational and/or electronic 
damping are included. 

With the above considerations, here we can treat (3 
as a fitting parameter, and approximate the higher rate 
of an adatom on the illuminated island edge by qh = 
v c exp (— kT ^C AE ) when crossing an island corner to a 
neighboring edge. The average frequency for an adatom 
on an island edge of length m to cross the island corner 
with rate q c can be well approximated [OJ by ^q c when 
V c > V e . Assume that the numbers of adatoms on the 
x + , x~ , y + , and y~ edges at time t are n+, n~ , n+, 
and riy , respectively. Because of the higher island cor- 
ner crossing rate for adatoms on the illuminated island 
edge, there must exist a net adatom current, ji, which 
transfers adatoms from the y + edge to its neighboring 
x + or x~ edge. As a consequence, one has n+ > n~ 
and n~ > n~ , which further induce another net adatom 
current, j 2 , carrying adatoms from the x + and x~ edges 
to the y~ edge. All together, the growth of the island is 
described by 
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Variations of the island width, Am^a, and island length, 



Am y a, are given by Am, 



and Am, 



The ratio r = -All = Am y /Am x predicts the evolution 
of the island shape. A square-shaped island will remain 
a square if r = 1 (stable growth), and will elongate per- 
pendicularly to the incident direction if r > 1 or parallel 
to the incident direction if r < 1 (unstable growth) . 
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FIG. 2. Temperature dependence of island growth on an 
fcc(100) surface with different contributions of condensation 



energy. 
m° = 2,m° y 



The initial island size is m. 



2 in (a) and 



= 8 in (b). Atoms are deposited at a grazing 
angle of a = 85°. Solid, dotted, dashed, and long dashed lines 
correspond to (3 = 0,0.0001,0.0006, and 0.05, respectively. 
Curves for (3 = 0.05 are truncated for r > 5 in (a) and r > 10 
in (b). 

Because AE reaches its maximum when a approaches 
90°, the effect of the condensation energy is most pro- 
nounced in the case of grazing incidence. In the following, 
we first focus our attention on the growth at grazing inci- 
dence. Typical behaviors of r versus T for different (3 are 
shown in Fig. 2, where the growth parameters are F = 
O.lMLs" 1 , v s = v c = 10 12 s-\ V s = 0.35eV, V e = O.leV, 
AV = 0.3eV, U = 5cV, and t = 0.1s. Fig. 2(a) is for 



an initial square island of size m x 



2. Fig. 2(b) 



is for an intermediate rectangular island with m x = 2 
and m y = 8. As expected, without contribution of the 
condensation energy, i.e., (3 = 0, the ratio r remains a 
constant r = 1 when T » T c and r < 1 when T < T c , 
where T c ss 190K is the freezing temperature for island 
corner crossing defined by tq c (T c ) ~ 1. For a non-zero (3, 
we find three different r — T behaviors upon changing the 
energy transfer coefficient (3. In the case of a very small 



j3, e.g., (3 = in Fig. 2, r increases with decreasing 

temperature until it reaches T c . Further decreasing T 
rapidly turns r to zero. A notable feature of r in the 
temperature regime T > T c is l<r<2. For larger (3, 
e.g., (3 = 6 x 10~ 4 in Fig. 2, when T approaches T c , r in- 
creases from r = 1 at high temperatures to a ratio r = 2. 
After that, a highly anisotropic growth with r > 2 ap- 
pears until the temperature reaches another critical value 
T a w 170K. When T further goes down, r drops to zero. 
We identify T a to be the freezing temperature for con- 
densation energy assisted island corner crossing defined 
by tqh(T a ) ~ 1. In the case of a much larger coefficient (3, 
e.g., (3 = 0.05, a steady growth regime characterized by 
r = 2 exists at moderate temperatures, and the strongly 
anisotropic growth mode r > 2 remains in the whole 
temperature range of T < T c . We note that the steady 
growth regime exists for all (3 > 0.05. Figure 2 also shows 
that the value of r is larger for a larger value of m y /m x 
(except for the steady state regime where r = 2), indi- 
cating that the elongation instability is more pronounced 
for already developed rectangular islands. 
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FIG. 3. Island growth on a Cu/Cu(001) surface at 
off-normal incidence and four different temperatures. 

The observed different behaviors of r for different (3 
can be rationalized based on the competition of the three 
rate processes: surface diffusion, thermally activated is- 
land corner crossing (without the assistance of the con- 
densation energy), and condensation energy assisted is- 
land corner crossing. In the temperature regime T > T c , 
frequent corner crossing with q c >> 1 and > q c in- 
duces a net current j\ transferring adatoms from the il- 
luminated edge to its neighboring edges. On the other 
hand, q c » 1 also results in a uniform adatom distri- 
bution on the x + , x~ , and y~ edges, with j 2 = 0. From 
Eqs. (2)-(3), we have r = 2/[l + (1 + 2m/n)£L where 
£ = (fy + Id - 2ji)/(2/ x + fy + ji). Because typically 
V s < V c , surface diffusion is activated when q c » 1, 
leading to f s » fd- Moreover, for a small (3, the net 
current is small with ji < (fd + /s)/2. Therefore, we 
have 1 < r < 2 for T > T c , as shown in Fig. 2 for 
(3 = 10~ 4 . By increasing (3, j\ reaches the maximum cur- 
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rent j\ = (f s +fd)/2 upon decreasing T toward T c . In this 
case, we find r = 2, as shown in Fig. 2 for (3 = 6 x 10~ 4 
and /3 = 0.05. In the temperature regime T < T c for 
small (3, corner crossing and surface diffusion are frozen 
with q c << 1, qh << 1 and q s << 1. The landing 
atoms from the direct flux fd accumulate on the illumi- 
nated edge, leading to w r < 1. For T a < T < T c , 
we have q c << 1 and q s << 1 but q^ » 1. Such a 
regime of strongly anisotropic island corner crossing leads 
to ji = (j s + jd)/2, which in turn results in strongly 
anisotropic growth with r > 2, as shown in Fig. 2 for 
(3 = 6 x 10~ 4 . Finally, for even larger /3, we have qu >> 1 
for any temperature below T c , leading to r > 2 for T < T c 
(see Fig. 2 for j3 = 0.05). The above discussions are valid 
as long as the islands are compact fl5| . 

Finally, through studying Cu/Cu(001) growth, we ex- 
amine how the aspect ratio r depends on the incident 
angle. Fig. 3 is obtained for various temperatures with 
the same flux F = 0.0042MLS" 1 as used in the experi- 
ments f| . The barriers are V s = 0.505 eV, V e = 0.265 eV, 
and AV = 0.29 eV, as suggested by experiments Q and 
embedded-atom model calculations |l5"| . The condensa- 
tion energy |j| is taken as Uq = 3eV and a reasonably 
small energy transfer coefficient, (3 — 8 x 10~ 5 , is used. 
Figure 3 indicates that, at room temperature, the elon- 
gation is not easily measurable for smaller a. Significant 
change of the island shape occurs at larger grazing an- 
gles (a > 80°). This explains why one usually observes 
square-shaped islands during deposition at normal inci- 
dence or small-angle incidence j4|j| but elongated islands 
at grazing incidence H. Furthermore, after integration 
of Eq. (2) to a coverage 9 = 0.5, we found that the as- 
pect ratio has a value r = = 1.05 for T = 250 K 
at a — 80°, as found in the experiments ]^). Moreover, 
Fig. 3 shows that r can be greatly enhanced by slightly 
decreasing the growth temperature, a prediction to be 
confirmed in future experiments. We also note that if (3 
increases weakly with a, the crossover shown in Fig. 3 
will be even sharper. 

In summary, through studying monolayer- high island 
shape evolution on an fcc(100) surface, we have shown 
that the condensation energy of deposited adatoms can 
play an important rule in controlling the island shape 
during epitaxial growth. In the case of deposition at off- 
normal incidence, the component of the deposition flux 
parallel to the surface provides additional atoms on the 
illuminated island edge and thus more condensation en- 
ergy. This leads to an enhancement of the mobility of 
the adatoms on the illuminated edge, and results in is- 
land elongation perpendicular to the incident direction. 
Such an island symmetry breaking phenomenon is most 
pronounced at grazing incidence. A comparison between 
the theoretical predictions and the experimental findings 
in Cu/Cu(001) growth confirms the general validity of 
the model. We have also found strong temperature de- 
pendence of the aspect ratio, and the existence of a well- 



defined incident angle above which the elongation insta- 
bility is most pronounced. These latter predictions are 
to be verified in future experiments. 
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